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Abstract 



Self-similarity of Burgers' equation with soine stochastic advection 
is studied. In self-similar variables a stationary solution is constructed 
which establishes the existence of a stochastically self-similar solution 
for the stochastic Burgers' equation. The analysis assumes that the 
i -Q . stochastic coefficient of advection is transformed to a white noise in 

C^ I the self-similar variables. Furthermore, by a diffusion approximation, 

the long time convergence to the self-similar solution is proved in the 
sense of distribution. 

^r^ . Keywords Self similarity; stochastic Burgers equation; diffusion approxi- 

^ i mation 
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2 : 1 Introduction 

(N 

The deterministic Burgers' partial differential equation for a field w{t, x) is 

K^- Wt = VWxx - WWx (1) 

. - -' and was proposed by Burgers ^ to help understand the statistical theory 

of turbulent fluid motion. Here w{t,x) is analogous to the velocity fleld and 
V represents the dissipative viscosity. To better model the randomness in- 
herent in the presumed chaos of turbulence, the following stochastic Burgers' 
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equation has been suggested [3 El [151 [13 [ISl [221 [21] and studied recently by 
many people [21[a[Sl[ini[ISl[2ni[2Il[2B]: 

wt = iywr,r,-wwx + h{t,x,w,Wa:) (2) 

where h(t,x,w,Wx) represents stochastic effects defined on a complete prob- 
ability space {fl,J-',F). On bounded domains the existence and unique- 
ness of global solution was studied by Da Prato et al. [20] when the noise 
term h(t,x,w,Wx) is white in time, and by Holden et al. [12] using a white 
noise calculus. On an unbounded domain the existence of a solution was stud- 
ied by a Cole-Hopf transformation by Bertini et al. [3] with h{t, x, w, Wx) an 
additive space-time white noise. 

We consider a family of solutions with special spatial-temporal form, 
namely the family of self-similar solutions, of the stochastic Burgers' equa- 
tion (|2|) on the unbounded real line with the particular stochastic advection 
h(t,x,w,Wx) = {wOx for some special space-time noise ((ty^) to be defined 
later. Here for the stochastic system, the self-similarity is in the sense of 
distribution which is defined later. The existence of self-similar solutions 
and the asymptotic emergence of self-similar solutions comprises the self- 
similarity of the stochastic Burgers' equation. Importantly for applications, 
the form of the stochastic advection, {w()x, is appropriate for globally con- 
served quantities w. Such stochastic advection is potentially of great interest 
in applications as it potentially illuminates some of the stochastic nature of 
chaotic turbulence in fluid flows. A thorough understanding of turbulence 
remains an outstanding challenge and researchers are increasingly invoking 
stochastic terms to model its effects in important environmental applica- 
tions [m [TH [391 6-g-]- We need to know how stochastic advection affects 
long term dynamics. 

Self-similarity is an important property of systems of physical interest, of 
which Burgers' equation is a special case. Many researchers have studies the 
existence of self-similar solutions of deterministic systems [5l [131 [211 [321 [13 
e.g.], and described the asymptotic behavior of self-similar solutions [21 [131 
[26l Uni e.g.]. But very little is known about self-similarity in stochastic 
spatio-temporal systems. 

We prove the existence and emergence of self-similar solutions, in the 
sense of distribution, for the stochastic Burgers' equation ([2|) for t > 1 . The 
stochastic advection h = {w()x for the stochastic Burgers' equation ([2|) trans- 
forms to a multiplicative white noise in the following self-similar variables. 
As in earlier research [21 [221 [23 e.g.] we introduce log-time and stretched 
space. 



and then define the stochastic field 

u{t, ^, u) = Vtw{t, X, w), tu E il . 

Straightforward algebra derives that in a useful sense, the spde ([2]) trans- 
forms in the similarity variables to 



du 



1, 1 



dT+{udW)^. (3) 



We call a solution w{t, x, u) to the stochastic Burgers' equation ([2]) a stochas- 
tically self-similar solution if the distribution of \/tw{t, x, oj) just depends on 
the self-similar variable ^ = x/\/t. By this definition, any stationary solu- 
tion u{^, u) to equation ([3]) is a stochastically self-similar solution of stochas- 
tic Burgers' equation ([2]). In order to construct a self-similar solution of 
stochastic Burgers' equation (|2]) we assume that W{t,^,u) is an L^(R) val- 
ued Wiener process defined on [Q, J-", P) with covariance operator Q which 
is detailed later. 

To construct a stationary solution of the transformed SPDE (jH]), we con- 
sider the system in a weighted space L'^{K) which is defined in the next 
section. First, by using energy estimates and the compact embedding re- 
sults of the weighted space, we show the tightness of solutions with initial 
value in the space L'^{K). Then the classical Bogolyubov-Krylov method [1] 
implies the existence of a stationary solution of the spde ([3]). Due to the 
multiplicative noise, the method showing the attraction of the stationary so- 
lution in the case of additive noise [S^ fails here due to the appearance of 
the unbounded term W. Instead we apply a diffusion approximation to this 
stationary solution. For the approximating model, which is a Burgers' type 
equation with a singular random perturbation (equation fHBl) ). attraction to 
the stationary solutions is derived by the method used for the case of additive 
noise. Then by the approximation, the attraction is passed to the stationary 
solution of ([3]). Here the most difficult part is to show the effectiveness of 
the approximation. We follow a martingale method to show the tightness of 
the family of stationary solutions of the approximating model. Then passing 
to the limit derives the attraction of stationary solutions of the stochastic 
advection Burgers' equation ([3]). 

2 Preliminary 

We consider the stochastic pde ([3]) in the self-similarity variables. For brevity 
we introduce the linear operator 



Denoting the weight function by K{^) = exp{,^^/4z/}, introduce the following 
weighted functional space for exponent p > 

and for positive integer exponent k 

H\K) = LeL\K):\\u\\l,^^^= J^ ||D"..||i.(^) < oo 

I 0<Q<fc 

Then the linear operator C is self-adjoint and generates an analytic semi- 
group S{t) on the space L'^{K) with the domain D{C) = H'^{K) ^. Fur- 
ther, the eigenvalues of the operator £ are 

k 
Afc = — - , A; = 0, 1,2, ... , 

with the corresponding eigenfunctions ek{C,), 

eo(0 = ^^exp{-eVM, e,(0 = Cfc%o(0, A; = 1,2,..., 

which forms a standard orthonormal basis of L'^{K) when we choose c^ as 
the constants such that ||efc||i;,2(^) = 1 . 

In the following we denote by Ec = span{ei(^)} and 

Es = E^ = Le L\K) : j^u{0 rf^ = o| . 

We also denote by Pg the linear projection from L'^{K) to Es. 

The following are some important basic properties on these weighted 
spaces [20] . 
Lemma 1. 1. The embedding H^{K) C L'^{K) is compact. 

2. There exists C > such that for any u G H^ {K) 

Jr Jm. 

3. For any u G H^{K), 

-^ Jr Jr 



4- For any u E Eg , 

{Cu,u) < --\\u\\l^K) 

5. Ifue H\K), then K^/^u e L° 



6. For any q > 2 and e > there exists constants C^^q > and R > , 
such that for any u e H\K) n L^^^(M) 

II l|2 ^ II ||2 I /-Y II ||2 

||M||i2(^) ^ ^\\U^\\l^(K) + '-^e,g||M||L9(B(o,_R))- 

Remark 1. By item [3] in the above lemma, in the space H^{K) we define the 
norm ||Vu||l2(x) which is equivalent to ||M||Hi(i<'). 

Further, by the spectral properties of the linear operator C, we define 
{—C + 1/2)''' for any 7 G M [38]. Then define the Sobolev space H'^{K), 
for any 7 e M, as V{{-C + 1/2)^/2), the domain of (-£ + 1/2)^/2^ By the 
embedding theorem [38], H'^'^{K) is compactly embedding into H'^'^{K) for 

7i > 72 • 

We make the following assumptions on the stochastic force. 
Assumption 2. 1. The stochastic force t\/trj = t\/t{wQx is written, in 
the self-similar variables, as {uW)^. Here W is an L'^{K) -valued Wiener 
process with covariance operator Q such that 



QV{0= qi^XMOKiOdC forany^eL\K), 

Jr 

with g(^,C) = Q'(C)0 positive, and 

[ [ q{^,C)K^{OKdC)d^dC<oo. 

Jr Jr 

The covariance Q shares the same eigenbasis as that of the operator C 

2. W^{t, ^) is an L'^{K) -valued Wiener process with covariance operator Q' 
such that 

QV(0= [q'i^XMOdC 



with q'iCC) = 9'(C)0 positive, and 

[ [ q'i^X)K{OKiC)d^dC<oo. 

Jr Jr 

Furthermore, 

TrQ < 00 and TrQ' < 00 , (4) 

andq{0:=q{^,O^H\K), 

||g||L°°(]RxK) and llg'llL-'(MxR) are small. (5) 



From the above assumptions, q'{^,C) = 9^9(^q{^,C) ^^^ the Wiener pro- 
cess W has the series representation 



fc=0 

where {Pk}k are independent standard Brownian motions. 
Remark 2. The special assumptions on r]{x,t,w) does not exclude the exis- 
tence of self-similar solutions for other cases. 

For later purposes we consider the spde ([3]) on the canonical probability 
space {Qq, J-", {JvjreR, P) which consists of the sample path of W{-,uj) in the 
space C(]R, L^(i^)) with Wiener measure P on Qq [1\ . To be more precise, 
W is the identity on fl, with 

no = {we C{R, L\K)) : w{0) = 0} . 

Let 9r : (f2o, J-", {-7v}reM)IP) ~^ (^O)-^? {-^rjreRjIP) be a metric dynamical 
system (driven system), that is, 

• 6*0 = id, 

• 6^6s = Or+s ioT all s, r G M , 

• the map (r, u) H- drOo is measurable and ^t-P = P for all r G M . 

On this canonical probability space Qq, we choose 6^ to be the Wiener shift [1] 

Orui-) = uj{- + t) -uj{t), teR, ueQo, (7) 

which preserves the Wiener measure P on Q^. Furthermore, 6r is ergodic 
under Wiener measure P. Writing W{t,C,) as W{t,^,u) to exphcitly show 
the dependence on a; G fio , then 

ly (-, e, 9rCo) = W{- + r, e, CO) - W{t, ^, u). 

Remark 3. Notice that Qq is different from Q. However, by the self similarity 
transformation, for any a; G f2 , there is a sample path of W{t, ^), then there 
exists an uq G Qq that represents the sample path of W{t,^) in Qq. The 
reverse is same. So in effect the space Qo is the same as Q. 

Recall that a random process {u{T)}r>o is said to be stationary if its joint 
probability distribution does not change when shifted in time r \i\. For the 
SPDE ([3]), to construct a stationary solution it is convenient to consider the 
transition semigroup associated to equation ([3]). We define {P^}^>o on the 



6 



space consisting of bounded continuous functions tp : L'^{K) fl L°°(]R) — ?► R 
by 121] 

(P.V^)(M°)=E^(M(r;M°)), (8) 

where u{t; m^) is the solution of equation (^ with initial value m° G L'^{K) fl 
L°°(]R). Denote by M. the space consisting all probability measures on the 
space L'^{K) fl L°°(]R) and endow Ai with the topology of weak convergence. 
Define the dual semigroup {P*}t>o acting on M. as 



L2(i<')nL°°(K) JL2(X)nL°=(R) 

for any ^ E Ai and bounded continuous function if) : L'^{K) fl L°°(R) — )■ 
M. If /^(m*^), the distribution of initial values m", equals /i, then P^/i is the 
distribution of the solution u{t;uo) jSHi Proposition 11.1]. Sometimes M. is 
too large, so we need the smaller space 

M2 = \ fJ, E M : / \\u\\L2fK)H{du) < 00 > . 

A probability space fi E Ai is said to be a stationary measure for the 
stochastic Burgers' equation ([3]) if 

P*fi = /i , for all t > . 

The following property of stationary measure is useful [291 Proposition 11.5]. 
Lemma 3. If fi E Ai is a stationary measure for ^ and the initial value vP 
is J-Q measurable with C{vP) = fi , then the solution process u{t; u^) is a 
stationary solution to the stochastic Burgers' equation ^. 



3 Existence of self-similar solutions 

By definition, a stationary solution to the spde ([3]) is a stochastically self- 
similar solution to the stochastic Burgers' equation ([2]). Next we construct 
a stationary solution to the spde ([3]) from any initial value Uq G L'^{K) fl 
L°°(M). 

For any r > , in the mild sense, the transformed stochastic Burgers' 
SPDE (ISD is written as 



u{t) = S{t)uq+ S{T-s)u{s)u^{s)ds+ S{t - s){u{s)dW{s))^ . (9) 
Jo Jo 

Then by the standard method for the existence of mild solutions to spdes [29] 
we obtain the following theorem. 



Theorem 4. Assume Assumption IB holds. For any T > and initial value 
uq G L'^{K) n L°°(]R), there is a unique solution u{t,^) to spde ^ in 
L^in,CiO,T;L\K))nL\0,T;H\K))). 

We construct a stationary solution by the Bogolyubov-Krylov method. 
For this we prove the tightness of the distribution of u{t, ^), {C{u{t, 0)}t>o , 
in the space L'^{K). We need some estimates in the spaces L°°(M) and L'^{K). 

3.1 L°°(M) estimates 

We follow the approach for a scalar convection-diffusion equation [10] which 
was also applied to characterise solutions to a stochastic Burgers' equation 
with additive noise [37] . 
We introduce 

+ fl, M > 0, Jl, M< 0, 

sgn(M)+ = \ and sgn(M) = <^ 

10, M<0; 10, M>0. 

Then for u G L'^iM) with Wu{t) e ^^(M), the integral ^^u^^<j){u) d^ < 
(> 0) for any nondecreasing (nonincreasing) cf) G C^(R). By a density 
discussion the integral j^u^^sgniu)^ d^ < (> 0). Moreover, the inte- 
grals f^uu^sgn{u)'^ dC, = and f^{C,u^ + u)sgn{u)'^ dC, = 0. Denote by 



u 



± 



sgn(M)=''M. Let m = ||'Uo||loo(ir). Then multiplying sgn(M — m)"*" and 
sgn(M^ — m)"*" on both sides of Q and integrating on M x [0, r] with r > , 
the integral 

/(M(r,O-m)+t/e<0. 

Jr 

Therefore, u{t,^) < m for any r > 0. Similarly u{t,C,) > —m for r > 0, 
Then ||ii(T)||L°°(]R) < m for all r > . 

3.2 Estimates in the space H^{K) 

We first give a uniform estimate in the space L'^{K). 

Let u{t, ^) = Uc{t, ^) -|- Ms(r, ^) with Uc G -Ec and Ug E Eg . Then 

(iwc = , 

dus = [Cug — Ps{uu^)] dr + Ps{udW)^ . 

So 

Mc(r,0 =Mc(0,O = («o,eo)eo(0 = / uo{^) d^eo{^) 



which is totally determined by the mass of tha initial value, namely M 
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Now applying Ito's formula to ||ms(t)||^2(;^), we obtain 



^^rll^^Wlli^W < ~lhs\\%i(K)- {uu^^u) 



dr 



+ i 



H\\Iq + \HIq' +{{uW)i,u). 



Here 



and 



are the norms defined on the Hilbert-Schmidt spaces 



C2{Q^^'^L^{K),L^{K)) and C2{Q'^/^L\K),L\K)) respectively [29]. Notic- 
ing that ||g||Loo(RxM) and ||g'||L°o(MxR) are small enough (jS]), we have some 
positive constant c such that 

,d 



UJT 



Integrating by parts yields 



Jr 



By property E] in Lemma [T], for any e,e' > and q > 2, there exist positive 



constants C^, Cs/^q and R such that 



[ufK^d^ 



f uiK'l\ufK'l^d^ 

Jr 



< 



1 



2 

< c 



1/2 



UufeKdi 

JR 

L'^(K)\W\\l'2{K)\\u\\l^ 



ufKdi 



1/2 



< 3eC||«dli^(A') + 3a 



L^K) 



, ^ II l|2 
+ '^e',q\\U\\Lql^B{0,R)) 



\U\ 



|M||loo( 



with some positive constant Ce',q.R- Then 

|(mM5,m)| < [eC + £'C^||M||2^(igJ \\U(.\\l2^^K)+CeCe>,q,R\M'L--{R)■ (10) 

Then for any e and e' > , there are positive constants that we still denote 
by Cg and C^'^q^R for some positive q and R such that 

1 II / \ II 2 

< -c||m.||hi(K) + c\W\\hHK) + [f^C* + e'Ce\\u\\\oo(^^)] \\ui\\] 

+ CeCe' ,q,R\\u\\\^^(^) + ((mW^)5,m) 

< [-c + eC + e'Ce\\u\\L^(^^:^'\ ||ms||hi(/^) 
+ [c + £:C + £:'Cs||m||loo(i;)] ||uc||ii(/^) 

+ C'eC'e'.g.ijIkllioom) + ((Miy)^,^) . 



Il2(x) 



Now choose e and e' > small enough and since Uc = Meo(0, ll'"c||ii'i(i<') < 
CM for some C > , then by the Gronwall lemma 

E||M(r)||i2(^)nLoo(K) < Ri , for all r > , 

with some positive random variable Ri. 
Given any ri > , in the mild sense 



u{t + Ti) = S{t)u{ti) + / S{t — a)u{Ti -\- a)u^{Ti + a) da 

Jo 

+ f Sir-a)itiin + a)dWin + a))^. 



JO 

Then taking H^{K) norm in the above equation, for any 1 > r > , by the 
properties of the semigroup S{t) there is some positive constant C, which 
also depends on q and q', such that 

E||M(r + ri)|||i(^) 



1 ^' 



< C(l + -=) E||n(ri; 



Il2(x) 



By the Gronwall lemma, there is a constant R2 > such that 

for any ri > and < r < 1 . Then we have uniform i7^(i^)-estimates 
for u'^{t) for all r > . By the compact embedding H^{K) C L'^{K) we 
establish the tightness of {C{u{T))}r>o, the distribution of {u{r)}r, in the 
space L?'{K). Then the classical Bogolyubov-Krylov method [1] yields the 
existence of one stationary measure denoted by yU. Let u{t, ^) be the solution 
with initial value distributing as fi, then u{t, ^) is a stationary solution to the 
transformed SPDE ([3]) with distribution fj, (Lemma[3]), and by the construction 
of the stationary solution, we have u G H^{K). 

Theorem 5. Assume Assumption\^ holds. For any initial Uq G L'^{K) fl 
L°°(M), there is a stationary solution, denoted by u, such that u G H^{K). 
Further, there is a sequence Tn, with r„ — )■ 00 as n —)■ 00 , such that 

u{Tn) converges in distribution to u in L?{K) 

as n ^ 00 . Here u{t) is the solution to the spde ^ with initial value uq. 
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Next we show that the above convergence holds for any sequence r„ with 
r„ — 7- oo as n — 7- oo . It is impractical to follow the approach used for the 
stochastic Burgers' equation with additive noise p7]. As Section [1] states, we 
consider in the next section a Burgers' equation with random fast fluctuating 
advection which is an approximation to the SPDE (|3]). 

4 Burgers' equation with random fast fluctu- 
ations: self-similar solution and stability 

We consider the following randomly fluctuating advection in a Burgers' type 
equation 

ul = Cu'-u'ul + ^^iu'f]')^, u'{0) = uo, (11) 

where f/^(r, ^) = f]{T/e, ^) in which fj{T, ^) is the stationary Ornstein-Uhlenbeck 
process solving 

dri = -r]dT + dW{T,^). (12) 

Then 

W{r, Ofis, C) = ig(e, C) exp (-^^) (13) 

and ?7'^(r) G H^{K) for any r > . By the assumption on W^{t,^) we have 
that the process fj^{T,^) = f/^(r/e, ,^) which solves 

dT]^ = -7]^ dt + dWi:{T, (14) 

and also f]^{T) G H^{K) for any r > with 

mi{r,0rjlis,0 = ig'(e,C)exp (-^^) • (15) 

In the following we sometimes write u explicitly to show the dependence 
on a; G Qq. By the definition of 6',-, the stationary process f]{T) = f]{T, u) can 
be written as fj{6rOj) := 77(0, OrUj). 

We study the dynamics of the random differential equation (rde) flTTj) for 
fixed e > . First, for any r > , in the mild sense the RDE flTTl) is written 
as 

u'{t) = S{t)uq+ I S{T-s)u'{s)u\{s)ds + ^ I S{T-s){u'{s)f]'{s))^ds. 
Jo V ^ Jo 

(16) 

Then by the classical theory for the deterministic Burgers' equation [35] this 

theorem follows. 
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Theorem 6. Assume Assumption [H holds. For any T > and initial Uq G 
L'^{K) n L°°(M) the rde 07]) /ias a unique mild solution 

u' G L^iVl, L\0, T; H\K)) n C(0, T; L\K))). 



Notice that the RDE flTT|) is an evolutionary equation with a stationary 
force. So the tightness of its solution immediately yields the existence of a 
stationary solution. We next construct a stationary solution for the RDE flTTl) . 
which is attractive for any e > . We follow the discussion in section [31 

First, by the same discussion as in section [XT] 

||m^(t)||l°°(ir) < ""^ for all r > , 

with m = \\uo\\l°°(r) ■ Next we give a uniform estimate in r in the space H^{K) 
from the estimate in space L'^{K) for u*^ with any fixed e > . We follow the 
same discussion as in section 13.21 

Let u''{t, ^) = Uc('7", + '"s('^' with Uc G Ec and Us E Eg . Then 

dul = 0, (17) 

rf< = [Cul-Ps{u'ul)]dT + j=^{uY)^. (18) 

So 

Wc(^, = Wc(0, = (%, eo)eo(0 = / %(0 d^eoiO 



which is totally determined by the mass of the initial value M := J^^ uo{C,) d^ . 
Now multiplying m*^ in the space L'^{K) on both sides of the RDE fITT]) we 
have 

Integrating by parts, 



{u'u\,u') = -\ f {u'fK^di. 

Jr 



Then by the same discussion as for fITOl) we have that for any 6,e' > and 
g > 2 , there exist positive constants Cs, R and Ce',q,R such that 

Moreover, for any e > there is some positive constant C^ such that 

12 



Then for any e and e' > , there are positive constants that we still denote 
by Cg and C^i^q^R for some positive q and R such that 

1 ^11 (= / \ II 2 
2j^\Ws{r)\\LHK) 

_L /^ /^ IU,<:||4 

~r ^e^£',q,R\\"' IIl°°(R) 

•-"ell'/ 



1 



-e||2 II e||2 j^ || e||2 



< 



-h + eC + e'CJu' 



2 1^..,^ ..£11^ iIloc(k) + -^ 



m 



e||2 



sWmiK) 



+ 



eC + £ Ce||-u''||ioo(]K) + -^ 



m 



e||2 



cll//l(i^) 



^^ n n iu,<:||4 _L 1 /^ ii,f^' 
+ ^e^£',qM\^ IIl°°(r) + ~n^£\\n 



^e||2 



\mK) 



\u 



e||2 



Now choose e and e' > small enough and since u\ = Meo(O) II'^cII-H'H^) — 
CM for some C > . Then by the Gronwall lemma 

\\u'{T,u))\\L2(K)nL--m < RliOrUj), for all r > 0, 

with some positive random variable R^ which is tempered by the properties 
of f]'^. 

Now given any ri > , in the mild sense 



u 



(r + ri) = S{t)u'{ti)+ S {t - a)u' {n + a)ul{n + a) da 

Jo 

+ 7i/ S{T-a){u'{n + a)f]'{n + a))^da. 



Then, taking the H^{K) norm in the above equation, by the growth of the 
semigroup S{t) for some positive constant C^, 



\U [t + Ti)\\Hi(^K) 



< c(l + ^ \\u\t^)\\lhk) 



c 



1 



+ 



c_ 



\/t — a 
1 



e-i--)/^Rl^e„uj)\\u'{n + (j)\\m(K) da 



e (^-'^)/2 I II -.^^^ ^ a)\\Lo.ij^)\\u\T^ + a)ll^i(^) 



\/t - a 

+ ll^1li/i(X)||'"1|L°°(K)] C?cr- 
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By the definition of r]'^ , tfie tempered property of Rl, and the Gronwall 
lemma for integral forms, there is a tempered random variable R2 such that 

\\u'{r + ri)||^i(;,) < CRl{9,co){l + j=)\\u^{t^)\\l,^k) + R\{e,uj) 

for any Ti > and < r < ri . Then we have uniform H^ (if )-estimates 
for M^(r) in r. By the compact embedding H^[K) C L'^{K) and the classical 
Bogolyubov-Krylov method [I], a stationary measure exists, denoted by /i^ 
Next we show the emergence of the stationary measure /x"^ by showing 
it is attractive. We follow the discussion for deterministic systems [lOl [21] 
which was also applied to the case of additive noise in a stochastic Burgers' 
equation [27|. The following lemma is a key step. 
Lemma 7. For any Ui,o, 1*2,0 ^ L'^i^) fl L°°(M) with 



Jr 



Let u\{t,^) and u^ir.S,) he the solutions to RDE ( fTI]) with initial value Ui^ 
and U2fl respectively. Then the function 



is strictly decreasing almost surely. 

Proof. Let U'^ij,^) = M^(r, ^) — U2{t,^), then it satisfies the following linear 
equation 

Notice that for any solution u"^ of the spde (11 II) . let f '^ = mI , then 



i^e 

5 



By the same discussion as in section 13.11 and the construction of r] 

(«i(r,0+«^(r,0-e-^r^^(r,0)? 

is bounded by a random constant for any r > . Then for almost all fixed u G 
Qq , the result follows by the discussion for the corresponding deterministic 
system [IDIEI]. □ 



Now we study the attracting property of any stationary measure. Notice 
that now m^ is not a Markov process, we consider process {u'^,ri'^) which is a 
Markov process on space {L'^{K) fl L°°(M)) x H^{K). 
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For this we introduce the space VJl consisting all probability measures on 
the space (^^(ir) nL~(R)) x H\K) and endow the space Tl with the topol- 
ogy of weak convergence. Define the continuous Markov semigroup {'P^}r>o 
associating with {u^,r]^) on DJl as 

V*MA) = F{{u^{t, ■),v'{r)) eA}, ueWl, 

for any Borel measurable set A C (L^(i^) nL°°(]R)) x H^{K). Here, the solu- 
tion to equations flTT]) - flT2|) . with initial value distributing as u, is (u''(r, ■), ri'^ij-)). 
Further, we introduce the following subspace of 9Jl; define 

m.2 = \^lem■ / (lhlli2(A') + MmiK))^^{d{u, //)) < oo I . 

L J{L'^(K)r\L^{M.))xm{K) ) 

Denote by v'^ the distribution of fj'^ on H^{K). Then by the construction 
of /x" G M.2, we have P*''(/2'', z/*") = (/i",//''), r > 0, which is a stationary 
measure of V^. Notice that r^'^, which is exponentially stable, is the unique 
stationary solution to equation flT2l) . So we just consider the projection 
of P;% which is denoted by P;% on L'^{K) n L°°(M). Then P*^Ji' = /i^ 
Moreover for any u = {fi, v) G OJI2 , there is a stationary measure u such that 
P*^u — >■ u as r — )■ 00 if and only if for any /i G A^2 , there is a /i" G A^2 such 
that P*''iJi — )■ /i^ as r — )■ 00 . 

We next show that for any /i G A^2 there is a unique stationary measure 
Ji" E M.2 such that P^*"/! converges weakly to /i^ as r — )■ 00 . 

Associated with the solution to the RDE (ITT1) we choose /x G A^2 which 
has the form 

fi = 6M*fJ's, (19) 

where (5m is some Dirac measure on Ec and fj,s is supported on E's. Then 
consider the limit of P*'^iJ, as r — )■ 00 . First by the uniform estimates in r 
for m" in H^{K) we have a measure /i% and subsequence r„ with r„ — )■ 00 , 
ra — )■ 00, such that 

P;>^/i% n^oo. (20) 

We next show p,'^ is unique by Lemma [71 Suppose /2''^ is another stationary 
measure such that for some r^ — )■ 00 , n — )■ 00 , 

P;f/i^/i'% n-^00. (21) 

Denote by 'u'^(r, ^) and u"^{t,^) the two solutions of RDE (TTTj) with initial 
value ^^(0 cind 'U^(^), distributed as /i^ and /i'^ respectively. Then 
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By Lemma [71 the function 

is almost surely strictly decreasing in r which contradicts the stationarity of 
u'^ and u"^. Hence we deduce the following theorem. 

Theorem 8. Assume Assumption {^ holds. For any initial uq G L'^{K) fl 
L°°(M), the solution u'^{t,C,) to RDE ( fTI]) . with initial value Uq, converges 
in distribution, as r —t- oo , to u'^ in the space L'^{K) which is the unique 
stationary solution to SPDE / flT]) with 



Jr Jr 






Remark 4. By the construction of u 



M'(r)|U-(R) <m, E\\u'{t)\\1^,j,.<C for all r > , (22) 



for some constant C > 0. 

We want to pass the above convergence property to the stochastic Burg- 
ers' equation Q, that is we want to pass to the limit e — )■ in uMn the 
space C{[0,oo),L'^{K)). We give some estimates uniform in e in the next 
section. 

5 Some a priori estimates on finite time in- 
tervals 

This section shows the tightness of {u''}o<e<i in the space C(0, T; L'^{K)) for 
any T > . 

From equation flTSl) . by the chain rule, 

lj^\\<ir)\\hiK) < -U^'ir)rH.^K) + {^'irH{r),u^{r)) (23) 

1 ^e/\„,e/\ I 1 -e. 



Define the two integrals 

Ilir) = -^ r{uUsms),u^{s))ds, 



mr) = l=J\u%s)vlis),u%s))ds 
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Now by the factorization method, for some < a; < 1 , 

-^ r iuUs)nHs),u\s))ds=^^^^ riT-sr-^YHs)ds 
V^Jo « Jo 

where 



Y'{s) = -^ / {s - r)-'^{ul{r)f]'{r),u'{r)) dr 
V ^ Jo 



s-r)''{{u'{r))lf]'{r))dr 
is-r) " / [u\r)Yf]\r)K^d^dr 



1 

1 



2v^Jo 



s — r] 



[u\r)Yf]l{r)Kdidr 



Then for any T > , there is some positive constant Ct such that 



sup |/i^(r)|' <Ct / \Yl{s)Yds + CT \Y^{s)Yds. 

0<T<T Jo Jo 

We first consider Y-^. By the L°^(]R) estimates on u"^, and the construction 
of fj", 



S PS 



< 



E\Y^'{s 

-^11 {s-r)-^{s-p)-" I {u%r,Oynr,OK^d^ 
^ Jo Jp Jr 

X [{u\pXW{pX)KcdCdrdp 

Jr 

p)-"- f fw{r;m\pX)K^K^didCdrdp 

JR. JR 



m 



Jp 



[s — r) [s 



and similarly 



E\Y,^{s)\' 

1 

e 



E/ J is-r)-"is-p)-'^ l_iu^ir,0)%ir,OKd^ 



X I {u'{p,C)fvlip,C)KdCdrdp 
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for some positive constant Ci t- Then 



E sup \II{t)\ <CtCi,t. 

0<r<T 

By the same discussion for Jg, a similar expectation holds: 

E sup \^{t)\<CtC2,t, 

0<T<T 

for some positive constant C2,t- Then by the same discussion for flTOj) and 
the Gronwall lemma, 

E sup ||w^(r)||i.(^)+E/ \\u^{s)\\l,^j,^ds<CT 

0<T<T Jo 

for some positive constant Ct- Notice that in the mild sense 

M^(r) = S{t)uq + / S{t — a)M''(cr)'U^(cr) da 
Jo 

+ ^ [ Sir-a){u'{a)r]'{a))^da. 
V e Jo 



(24) 



Then for any T > t > 6 > , 

< \\{S{r)-S{6))uo\\mK) + 



(25) 



S{T-a)u'{a)ul{a)da 



LHK) 



1 

7^ 



S{t - a){u\a)f]'{a))i.da 



LHK) 



f [S{t -a)- S{5 - a)]u'{a)ul{a) da 
Jo 



L\K) 



1 

7^ 



[S{t -a)- S{5 - a)]{u'{a)f]\a))i. da 



LHK) 



By the L°°(]R) estimate of u^, and by estimate flMj) . the expectation 



E 



S{t — a)u^{a)u''^{a)da 



LHK) 



< E / \\S{r-a)u'{a)ul{a)\\L2^K)da 



< E / \\u^{a)ul{a)\\L2(^K)da 



< mE / \\uj{a)\\L2(K)da < mCTVT^-S . 
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LHK) 



Expanding by {ek]k and by ([I3D, 

S{T-a){u'{a)f]'{a))^d(y 

X e-^*("-^) [ u'{\,C)f]'{X,C){ekK)^dCdadX 

Jr 

I, J& Js JrJr 



2 

m 



' k 



X {ekK)^{ekK)^ dC, dC da dX 

for some positive constant Ct- By the strong continuity of the semigroup S{t) 
and a similar discussion, from (125|) the expectation 

E\\u\T)-u'{5)\\L2iK)<CTV^^. (26) 

Now we need the following lemma [25j. Suppose Xi and X2 are two 
Banach spaces. Let T>0,l<p<oo, and i3 be a compact operator 
from Xi to X2] that is, B maps bounded sets of Xi to relatively compact 
subsets of X2. 

Lemma 9 ([2S])- -^^^ H be a bounded subset of L^{0,T] Xi) such that G = 
BH is a subset of L^(0, T; X2) bounded in L^{0, T; X2) with r > 1 . // 

lim ||n(- + a) — u(-)||lp(o,t;A'2) = uniformly for u E G , 

then G is relatively compact in Lp{0, T; X2) (and in C(0, T; X2) if p = +00 j. 
Let Xi = H^{K), X2 = L^{K) and B be the embedding from Xi to X2 , 
by Lemma [9], from estimates f l2^ and f l26|) we obtain the main theorem of 
this section 

Theorem 10. Assume Assumptions^ holds. For any T > , the distribution 
o/{u''}o<e<i is tight in the space C(0,T; L^(K)). 

6 Diffusion approximation 

This section determines the limit of m^, the solutions of the RDE (ITT]) , as 
e — )■ . We first show the tightness of u^ in the space C([0, 00), L'^{K)) by 
Theorem [TOl Then we determine the limit of u"^ in the space C([0, 00), L'^[K)) 
by a martingale approach. 
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6.1 Tightness in space C{[0,oo),L'^{K)) 

We need the following result on the tightness of a family processes fTR The- 
orem 3.9.1]. 

Lemma 11. Let X be a Polish space and let {X'^}o<e<i be a family of pro- 
cesses with sample paths in the space C{[0,oo),X). Suppose that for any 
6 > and T > there exists a compact set Tg^T C X such that for all 
0<e< 1 

F{X'{t) e Ts,T forO<t<T}>l-6. (27) 

Then {X"^} is tight in the space C([0, oo), X) if and only if {F{X'^)}o<ce<i is 
tight in the space C([0,oo),M) for any F G Cb{X) the space consisting all 
continuous and bounded functions on X. 

Remark 5. We do not need to verify the tightness of {F(X^)}o<e<i for all 
F e Cb(A'). One just needs to verify the tightness for all F in a dense 
subset of Ch{X^ in the topology of uniform convergence on compact sets [TT1. 
Theorem 3.9.1]. 

By Theorem [TUl the pre-condition flTTI) in Lemma [TT] holds. Next we show 
the tightness of {F(u"')}o<e<i in the space C([0, oo)) for any F G Ci,{l?{Ky). 
We follow a martingale approach. Continue to let A" be a Polish space and 
{X'^}o<e<i be a family of processes valued in the space C([0, oo), X^ adapted 
to the filtration T\. Let £^ be the Banach space of real valued J-'^-progressive 
processes with norm ||1^||£e = sup^>o^E|^(''")l- Let 

M" =[ (y, Z) e £' X £^ Y{r) - f Z{s) ds is J'^-martingale l . (28) 

Then the following lemma applies [HI Theorem 3.9.4]. 

Lemma 12. For any bounded continuous function F on X with bounded 

support, and for any S > and T > 0, there is {V, Z'') G A^*" such that 



< S (29) 



and 

supE [||Z'']|ip(o,T)] < oo /or some p G (1, oo]. (30) 

Then {F(X^)}o<e<i ^s tight m C([0, oo),M). 

By Remark [5] we just need to show the tightness of the following family 
of real valued processes |22j : 



lim sup E 

e-S>0 


sup Y^{t) - 

Te[o,T] 


- HX^r)) 



{/((«'. <^))W<i 
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for any </? G V{R) and twice differentiable compactly supported functions /. 
From the RDE fITT]) 



fi{u^iT),^))-fi{Uo,^)) 



(31) 



f'i{u^is),^)){u^{s),C^)ds- / f'i{u^is),^)){u^is)uUs),^)ds 



+ ^ / /'((n^(.),v.))((n^(.)r/^(.))5,^)cis. 
v^ Jo 

One can see that the singular term in the above equation is difficult to treat. 
To treat this term we follow a perturbation approach developed by Kush- 
ner [22]. Let J-"^ be the a-algebra generated by {^"^(s) : < s < r}. Then 
introduce the process 



F^ir) 



=E 



f'i{u^iT),^)){iu\T)r^^is))^,^)ds 



Tt 



(32) 



For the process -Ff (r) the following lemma holds. 
Lemma 13. Assume Assumptionl^ holds. Then 



FKr) = V~ef{{u'{T),vmu^{Tmr)h^)- 



(33) 



Furthermore, 



for some positive constant C, and for any T > 

E sup \FI{t)\ ^0, ase^O. 

0<T<T 



Proof. The equality fl33|) is implied by flT3|l and the property of conditional 
expectation. Then by the L°° bound on u*^ and the estimates on f/*^, 

E|Fi^(r)| < V~e\\f'\\L^iR)h'{r)\\L^iU)nmr)\\LHK)y^ + I^^l^^k) 
< C^t 

for some positive constant C. Further, by the maximal estimate on stochastic 
integral [2]?| Lemma 7.2], for any T > 

E sup U{T)fj^.(K.<CT 



0<T<T 

for some positive constant Ct- Then by 



E sup \F^iT)\ ^0, as e^O. 

0<T<T 



The proof is complete. 



D 
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To apply Lemma [T^ we first construct [V-, Z"^) G 2!^. For this we intro- 
duce the operator A'' defined by 

A7(r) = P - lim i [E/(r + 5) - /{r) \ T^] (34) 

for any J-"^ measurable function / with sup^E|/(r)| < oo . Then Ethier and 
Kurtz's proposition fHX Proposition 2.7.6] yields that 



fir)- / AJis)ds 
Jo 

is a martingale with respect to J-"^. Now define [V^, Z"^) as 

Y^ir) = fiiu^ir), v)) - Flir), Z^{t) = A^Y^{r). 

Then we establish the following lemma. 
Lemma 14. 

Z^r) = /'((«^(r),^))(n^(r),£^) + /'((^^(r),^))(i(n^(r))|,v.) 
+ fiiu^r), ^)){{u%T)nr)h iVi + l^ViWir)) 
- Vefiiu^ir), ^)) [{u^r), C{{^^ + \i^W{r)) 

Proof. By (jSI]), 

A^f{{u^{r),^) = f'{{u^{T),^)){u^{T),Cv) 

-f\{u\T),v)){u\r)u\{T),v) 

Now consider A'^F^. By (l33l) and the construction of f^*^ , 

nFi{T + 5) \ :fj,] 

= ^emf\{u\T + 6\^))-f\{u\T),^))] 

x(n^(r + 5)r(r + 5),0e + Kv5)l^;} 

+ v^/((«^(r), ip))n{[u\r + 5)- u\tM{t + 5),ip^ + \iip) \ :F^] 

+ v^/((n^(r),^))(«^(r)r(r)e-^/^^5 + \iv) . 
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+ —fi{u^iT),^Mu^irmr)hv) 



Then 

1 

This completes the proof. D 

Now by the above construction of {V^, Z^), 

Then by Lemma fT3| 

hmE sup \FI{t)\ =0. 

"^^0 0<r<T 

Furthermore, by the L°°(M) estimates on u'^{t), 

supE||Z^(r)|U2(o,r)<oo. (35) 

By Lemma [121 we thus deduce the following theorem. 

Theorem 15. Assume Assumption\^holds. The family of processes {M'^}o<e<i 

is tight in the space C([0, oo), L'^[K)). 

Remark 6. By (1221) . all discussions in this section and section O hold for u*^. 

So the above tightness result also holds for {'u''}o<e<i . 

6.2 Pass to the limit e — > 

We show the convergence of m^ as e — )■ and determine the limit. For this 
we introduce a limit martingale problem and show any accumulation point 
of {u^} is a solution to this martingale problem. By the convergence result 
of Walsh [36. Theorem 6.15], we just need to consider finite dimensional 
distributions of {(uS v^i), . . . , (u% V'n)} for any </?!,..., yj^ G I^(M); that is, 
we just pass to the limit e — )■ in the following equality 



E 



Y'{t) -YHs)- I Z'{r)dr 

s 



h{{u^{ri), (fi),..., (M'(r„), yp„)) > = 

(36) 
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for any bounded continuous function h and < ri < ■ • ■ < r^ < T with 
any T > . Denote by u one limit point in the sense of distribution of m'^ 
as e — )■ in space -D([0, oo), L'^{K)). By the estimates in Lemma [13] and the 
construction of V, the hmit of Y''{t) — Y''{s) is f{{u{T), ip)) — f{{u{s), ip)). 
Now we treat the integral term. First denote by Zl{-), A; = 1,2,3,4, the 
first four terms of Z'^{-) and by Z^{-) the last two terms in Z'^{-). Then, in 
distribution as e — )■ , 



J s 

and by the estimates on -u^ in section \5\ and estimates on f/^ in section [21 

E / \Z'^{r)\dr ^ . 

J s 

Notice that {u^f is bounded in ^^(0, T; L'^{K)) for any T > and by the 
tightness u*^ in the space C(0, T; L'^{K)), m^ converges almost everywhere to u 
on [0, T] X R, then by the L°°(R) bound on m^ and u, we have in distribution 
for any T > 

((m')^ ¥?) ^ (m^ V?) ase^O. (37) 

So in distribution as e — )■ 

j Zl{r)dr ^ I f\{u{r),ip)){u{r)u^{r),ip)dr. 

J s J s 

Next we treat terms Zg and Z^. For any u G L'^{K) define the bilinear 
operator S(m) such that 



(S(«)^,^)= / / «(0«(C)g(e,C)(^(0^(e)Mv5(C)i^(C))crferfC (38) 

Jr Jr 
for any ip G I?(M), and define the linear operator 

{A{u),^) = I fuiOq{^,0{^{OK{Ohd^ (39) 

Jr 

+ -, fuiOq{^,Oi^iOK{Ohd^. 



For this operator S the following lemma holds. 

Lemma 16. For any u e H^{K), let B{t,^) = {u{^)W{t,^))^ , then B is 

an L'^{K) valued Wiener process with the covariation operator S(m). 
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Proof. The proof is direct. By (Q, 

k=l 

Then by the representation of q{^, (), 

= EU(0$^v^e,(0/3fc(r)) ( «(C) $^ v^Q(C)A(r) 

\ k J ^\ I / ^ 

+ u{0u{0qui^,0. 

By the definition of S(m) this proves the lemma. D 

To pass to the hmit e — )■ in Zg and Zl we apply again the perturbation 
method [22]. Let 



F^r) 



1 t,^\2 



((«^(r),^)) / E {u^{r)ns),n + l^^) 



-(S(n^(r)V,^)|j-; 



ds 



and 



-{Aiu^{T)),^)\r^ 

By the construction of f]^ and fjj , 

-2(s-r)/e-e 



F4(r) = n{u\T),v)) E {{u\T)n\s)%{v^ + \^^W{s)) 



ds . 



-2{s-r)/e^ 



(40) 

(41) 
Then 

i^sM = f r ((«1r), v^)) [(«^(r)r (r), y.^ + l^^)' - (S(«^(r))v,, ^)] (42) 

and 

Flir) = f /'((«1r), (p)) [{{u^{T)nr)h in + 1^^^)^ M) - {A{u^{r)), v)] . 

(43) 
Then by the estimates on u'^{t), f]'^{T) and i]Ut), direct computation yields 
this lemma. 
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Lemma 17. As e ^ , 

supEF|(r) = C(e), and supEF|(r) = 0{e). 

T>0 T>0 

Now following the same discussion as in Lemma [H] and f HU]) -f HT]) . we 
have the following lemma. 
Lemma 18. 

A'Fiir) = r((«^(r),v.)[(S(«^(r))^,y,) 
and 



with 



as e — 7- . 



supE|i?^(r)| = C(e) and supE|i?^(r)| = 0(e) 

T>0 r>0 



Proof. This is similar to the discussion in the proof of Lemma [TH First we 
have for any 6 > 

E[F|(r + 5) I J-;] 
= fE[(r((«^(r + 5),^))-r((^^(r),^))) 

+ f r ((«^(r), y,))E [((«^(r + 5) - w^(r))r (r + 5), y.)^ 

+ U"{{u'{r), ^))E f(n^(r)r^^(r + 6), ipf - (S(«^(r))^, ^) | ^^ 
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and 



E[Fl{T + 6)\r^] 



m 



X {{{u^ir + 5)n\T + 5))^, {^^ + \ivW{T + 5)) 
-{A{u^{T + 5)),v))\r^_ 

+ f /'((«^(r), vm[{\{u\T + 5)- u^irWir + 6)]^, 

(cp^ + ie^^)r (r + 6))- {Aiu^ir + 6)) - A(«^(r)), y,) | J"; 

+ f /'((«^(r), v9))E[((«^(r)r(r + <5))^, (y.^ + ie¥')r(r + 5)) 

Then by the definition of A"" and (l40l) - fHTl) . direct computation yields the 
result. The proof is complete. D 

Now we have the following J-"^ martingale 

M'{t) = f{{u%r),^))-Fl{T)-F^{T)-F:iT) 

f'iiu^ir), ^)) \{u\r), Cif) + {\{u^{r)f, ^^ + \i^) 

dr 



+ {A{u\r)),v) 
f"{{u^{r), ^))(E(n^(r))^, ^)dr + R%t) 



where 



mr) 



[Zl{s) + Rl{s) + Rl{s)]ds 



with E|i?^(r)| = C(e) as e — )■ 0. Now passing to the limit e — )■ , the 
distribution of the limit u solves the martingale problem 



M{r) = f{{u{T),^))~ / f{{u{T),^)) 



(44) 
(M(r), C^) + (|(n(r))2, ^^ + l^^) + {A{u{r)), v^)] dr 
f"{{u{r), ip)){i:{u{r))ip, ip)dr 

which, by Lemma fT6| is equivalent to the martingale solution to the spde [27] 

Ur = Cu — uu^ + \{uq)i^^ — \{uq)^ + {uW)^ (45) 
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for some new Wiener process W with the same distribution as that of W. By 
the general theory of SPDEs [21] , the martingale solution to the spde f H5]) is 
unique in the space L'^{0, T; H\K)) n C(0, T; L^{K)) for any T > . Then 
we deduce the following theorem. 

Theorem 19. Assume Assumption\^ holds. The solution of rde / flT]) . u"^, 
converges in distribution in the space C([0, oo), L'^{K)) to u which solves the 

SPDE (P^. 



Notice that the SPDE f H5|l is different from the stochastic Burgers' equa- 
tion (|3]). For this we consider the following rde 

< = Cu' - u'u\ - \{u^q\^ + i(MV)€ + T^{'^'n% ■ (46) 

By the assumption that g(^) G H'^{K) and ([5]), the extra terms \{u''q)^^ and 

\{u''q')^ do not change the estimates in sections [SHBl Then we derive this 

theorem. 

Theorem 20. Assume Assumption [H holds. The solution of the RDE ^JW 



converges in distribution in space C([0, oo), L'^{K)) to a solution of the SPDE ^. 

Assume u*^ is a solution of fH6|) with initial value uq G L'^{K) fl L°°(R) 
and u*^ is unique the stationary solution such that 

u'^ij) converges in distribution to u'^ij) in L'^{K) as r — ;■ oo . 

Now assume u is the limit of m^ in distribution. Then for the solution u of 
the SPDE ([3]) with initial value Uq 

u{t) converges in distribution to u{t) in L'^{K) as r — ;■ oo . (47) 

Then we draw the following main theorem. 

Theorem 21. Assume Assumptionl^ holds. For any solution of the stochas- 
tic Burgers' equation ^ with w{l,x) G L'^{K) fl L°°(]R), there is a unique 
self-similar solution w{t, x) such that 

vtw{t, x) — Vtw{t, x) — )■ , as t —> oo in distribution in L^iWj. 
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